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+3 , Abstract. Upcoming Weak Lensing (WL) surveys can be used to constrain Dark 

Energy (DE) properties, namely if tomographic techniques are used to improve their 



sensitivity. In this work, we use a Fisher matrix technique to compare the power of 
CMB anisotropy and polarization data with tomographic WL data, in constraining 
DE parameters. Adding WL data to available CMB data improves the detection of 
all cosmological parameters, but the impact is really strong when DE-DM coupling 
is considered, as WL tomography can then succeed to reduce the errors on some 
parameters by factors > 10 . 
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1. Introduction 

> : 

^ . The first data system requiring Dark Energy (DE) concerned cosmic acceleration, 
detected through high-redshift supernovae [JJ. CMB [2] and deep sample [3] data 
supported also the DE case, showing that the density parameter for non-relativistic 
matter f2o,m ~ 0.3, while the total density parameter Qq ~ 1 . 

In the most popular scenario, DE is ascribed to a cosmological constant A. 
Alternative options include a self-interacting scalar field, <fi (quintessence or dynamical 
DE [HE]) and modifications of General Relativity [5J. 

It is known that models with A (ACDM) apparently accommodate all available 
data systems. The problem is the physical origin of A, which can be a false vacuum; 
this however causes well known fine tuning and coincidence problems. 

The former problem is partially eased in dynamical DE (dDE) scenarios, when self 
interaction is due to a tracking potential V((f>) [7j. If V((f)) is SUGRA [5] , the fit with 
data is at least as good as for ACDM [8]. 

In the attempt to ease the coincidence problem, DM-DE interaction (e.g., (9j [TO] ) 
was also considered, yielding an energy transfer between the dark components, so 
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allowing a (quasi)-parallel scaling of DM and DE from a fairly high redshift until the 
present. While laboratory data set no significant constraint on DM-DE interaction 
strength, parametrized by ft (see below), recent works placed constraints on possible 
couplings, by using SNIa data [11] or the redshift evolution of the Hubble parameter, 
H [12]. Accordingly, ft > 0.12-0.15 [131 [H] is hardly consistent with observations. 

Unfortunately, such a low coupling level no longer eases the coincidence 
problem [IB] , but, once the genie has come out from the lamp, it is hard to put it 
back inside. The point is whether low values of ft, as allowed by current data, can 
interfere with future data analysis. In particular, when we allow for non-zero ft, how 
do errors on other parameters behave? 

In this work we tried to answer this question by using a Fisher matrix technique. 
We considered two different models, set by similar values of cosmological parameters, 
without and with coupling. In the latter case, we took ft = 0.1 . Starting from these 
models, we evaluated the expected errors on cosmological parameters, as obtained when 
data concern just CMB anisotropy and polarization or include tomographic weak lensing 
(WL). 

As a matter of fact, in coupled models, the time evolution of the dark components 
is non-standard. If such models are considered in a Newtonian approximation, it is as 
though DM particles had a ^-dependent mass. Also for quite low ft's, this anomalous 
scaling leaves an imprint on both the expansion history of the Universe, and the growth 
of (matter) fluctuations, at the linear and non-linear levels (e.g. [16J). 

However, any detected evolution of H can be reproduced through a suitable redshift 
dependence of DE density pd e and state parameter w^, when <fi approaches m p (the 
Planck mass). A risk is that, if matter and dark energy are coupled, fitting observations 
leads to an estimate of a phantom equation of state (w^ < — 1), even if > — 1 at all 
redshifts [T7] . 

In principle, this risk can be excluded if the redshift dependence of the growth factor 
G(z) is also tested, through the increase in number and concentration of bound systems. 
Data providing information both on H(z) and G(z) are therefore able to discriminate 
between coupled and uncoupled models. Experiments, or combinations of experiments, 
probing H(z) and G(z) are then needed. 

CMB data, used to constrain coupling [131 [IB] , place only upper limits on ft. The 
analysis of Ly-a and the matter power spectrum of the 2dF and SDSS surveys [H] does 
not lead to great improvements. At the available sensitivity level, such data systems 
provide just weighted integrals of H(z) and G(z), which remain consistent with a rather 
wide set of options. 

On the contrary, gravitational lensing, alone or in combination with CMB data, 
was already shown to be a powerful tool for the analysis of DE. WL tomography probes 
the power spectrum P(k) at different redshifts and is thus well suited to constrain G(z) . 

In this work we aim to put these conceptual points on a more quantitative basis and 
to deepen the case of coupling, by performing a Fisher analysis of future WL surveys 
and CMB experiments. 
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The outline of this work is as follow. In Sec. [2] we review the basic properties 
and definitions of dDE models and WL, in Sec. [3] we show the results of the Fisher 
analysis, in Sec. H] we discuss them and in Sec. [5] we summarize our findings and draw 
our conclusions. 



2. Models and definitions 

2.1. Interacting Dark Energy 

We consider a cosmological model where the DE field <j) interacts with the cold DM 
component. The model requires the specification of the potential V((f>) and the function 
f((f>) characterizing the coupling. The equation of motion for then reads 

+ 3H<j)= -V e / f with V eff = V + p c . (1) 

Here dots denote ordinary time differentiation, H(a) = a/a and p c is DM energy density. 
In turn, its evolution is governed by 

p c + (3H + C<P) Pc = , with C{et>) = dl °^ f) . (2) 

This equation can be integrated and gives: 

pc(a) = PcocrVO) • (3) 

For / = 1 eqs. ([T]), ([3]) return ordinary dDE equations. The equations for the other 
components remain unchanged. In a generic coupled model, then, the ratio between the 
energy densities of cold DM and baryons is not fixed, but evolves in time according to 

However, it is always possible to define an effective DE component of density 

p e J/(a)=p c a~ 3 if(<i ) )-l]+p de (a) . (4) 

In general, p e J/{a) is not guaranteed to be positive and detecting p e JJ [a) < would be 
a clear indication that our description of the dark sector is not adequate. Lacking such 
clear giveaway, however, experiments probing H(z) can hardly discriminate between 
DE-DM interaction and an ad hoc DE component. If its effective state parameter 



w 



eff 



p e ff/p e ff^ i s < — 1 ) data may appear to require phantom DE (see [19] for 
discussion) . 

Coupling affects also fluctuation growth. In the Newtonian limit, i.e. well below 
the horizon, and neglecting the contribution of radiation, baryons and DM fluctuations 
grow according to the equations (e.g. [161 120]): 

5 h + 2H6 h = 4vrG(p b 5 b + pA) , . 

S c + 2HS C = 4ttG [ Ph 5 h + (1 + lP 2 ) Pc S c ] , 1 } 

where we defined (3 so that 
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Figure 1. The growth suppression rate g(a) = 5(a) /a for CDM (solid line) and 
baryons (dashed line) in a coupled model. Curves refer to a SUGRA potential with 
exponential coupling. It can be noticed that the cold dark matter evolution rapidly 
diverges from the standard cold dark matter solution, g(a) = 1 (dotted line), even 
well into the matter dominated era. For comparison, we also plot g(a) for the same 
SUGRA model with coupling turned off (dot-dashed line). 

while m v = G? -1 / 2 . Therefore, baryons and DM perturbations grow at different rates 
and, even soon after recombination, a growing mode 5 oc a no longer exists (see Fig. [[]), 
leading to a bias between baryon and DM perturbations. Analytical models of spherical 
collapse have shown that this differential growth results in a baryon-DM segregation, 
with baryons occupying the outer regions of collapsed objects [2T] . 

Eqs (JSJ) also show that the growth equations explicitly depend on both H and p c ; 
therefore, if we measure the growth from data, a possible anomalous scaling can no 
longer be masked though an ad-hoc definition of an effective DE density It is then licit 
to conclude that experiments probing the rate of growth of fluctuations are in principle 
well suited to test coupling between the dark components. 

In this work we are interested to combining WL and CMB data, which cannot be 
accurately described using current parametrisations [22J. Therefore, we follow a more 
conventional approach and choose the functional forms 

/<*>-p(V^)m 

The SUGRA [5] potential V((f>) depends on the slope a and the energy scale A. Fixing 
DE density today and A (a), however, determines a unique value of a (A). CMB, SNIa 
and deep sample data yield A < 10 3 GeV JS], in the absence of coupling. Here we focus 
on the SUGRA potential as it naturally arises in the context of Supergravity Theories 
and is an example of tracking potential characterized by a rapid time variation of the 
equation of state, when DE becomes dominant. Therefore, assuming a constant w for 
this class of potential may lead to misleading results. Different choices for the potential 
are clearly possible. 
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The coupling function /(</>) depends on f3, and 0o is the field value today. In this 
work we assume a constant j3 > (see however [19] for a different approach); data place 
the upper limit (3 < 0.12 — 0.15 [131 Ej- F° r reasonable values of the cosmological 
parameters and of A, we expect coupling effects not to be relevant for (3 < 0.01, so that 
the dynamically interesting values for the coupling lies in the range 0.01 < (3 < 0.10. 



2.2. Weak lensing 

Among the cosmological probes allowing the analysis of the nature of DE, the 
cosmological WL has been earning a fundamental role (see [231 EH EHl [26] for a thorough 
review). In fact, next generation WL surveys will cover a significant fraction of the 
sky and observe galaxies at deeper redshift, making WL a powerful tool to study the 
properties of the Universe. Furthermore, WL tomography will allow to significantly 
increase the cosmological information that can be recovered from such surveys (see [27] 
and [281 EHl [30] for details about lensing tomography). 

The power spectrum for the WL convergence between the zth- and jth- redshift 
bin is given by: 

p <«> (f > = (t)7 ^««* {rib)' z ) ■ (8) 

Here P^i,(£/r(0, z), z) is the non-linear matter power spectrum, at the redshift z and 
wave number k = £/r(0, z); the quantity r(z, z s ) is the radial comoving distance between 
z and z s , 

f Za cdz' 

r(«,«.) = jf j^, (9) 

while 

W i (z) = ^Q m0 F i (z)(l + z) (10) 
is the window function which weights different redshift bins according to the factor 

being Zi the i-th redshift bin. The function Ui(z) is the normalized redshift distribution 
of the source galaxies falling in the i-th photometric redshift bin, per unit solid angle, 

-i 



Ui(z) = Di(z) 



dz'DAz') 



(12) 



where 



(i+l) 
ph 



Di(z) = / dz ph n(z)p(z ph \z). (13) 

ph 

The last expression codifies the effect of errors in photometric redshifts of source galaxies 
[3~T] . In fact, due to the large number of observed galaxies with future surveys, one needs 
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to refer to their photometric redshifts, even if the determination of these redshifts may 
be not so accurate as the spectroscopic ones. In the model given by (fT3l) . the mapping 
between the photometric z p h and the spectroscopic redshift z is obtained convolving 
the overall galaxy distribution per unit solid angle, n(z), with a probability distribution 
p{z p h\z) in z p h at a given z. We choose a Gaussian function at each redshift for the 
distribution of photometric redshifts, i.e. 

(z — Zph)' 



p(z P h\z) = j= exp 



2.2 (14) 

while the overall distribution of source galaxies is chosen to have the parametrized form: 

where A, B and Zq are the parameters. One can easily check that equation (fT3l) turns 
into: 

Di(z) = -n{z) [erf(x m ) - erf(^)] , (16) 

with Xi = (Zpl — z)/y/2a z and erf(x) the error function. 

The cosmic shear power spectrum will receive a shot-noise contribution from the 
random intrinsic ellipticities of source galaxies and measurement error |32j. Thus, the 
observed power spectrum between redshift bins i and j can be expressed as: 

P$)(£) = P(^) + tS (17) 

• H 

where a t is the rms shear due to intrinsic ellipticity and measurement noise (we assume 
a e ~ 0.22 [271) and 

-2 



n g 



lamin 



1 7T 

60 180 



hi (18) 



is the average number density of galaxies per steradians in the i-th redshift bin, n g being 
the number of galaxies per square arcminute and hi the fraction of sources belonging to 
the bin. 

The covariance between the power spectra P°^(£) and P?^ n -\{^') is approximately 
given by 

cov [p^(£), p^ n) {i')] = {u - 5 ^ Mf [Pfa)W PffoW + PfaW p^W)] ( 19 ) 

where f s k y is the sky fraction covered by the survey and A£ is the bin width centred at 
i. The above expression assumes that the power spectrum in each multipole bin is very 
flat, in order to replace the value of the spectrum evaluated at the bin center with the 
average of spectrum over each bin (see Appendix A for the complete expression). In 
addition, we have not included the non-Gaussian term, due to the contribution of the 
shear trispectrum 
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2.3. Fisher's formalism 



The Fisher matrix formalism [35], [36], |37j provides lower limits to the error bars of the 
cosmological parameters one wishes to measure. The basic tool in Fisher's method is 
the likelihood function, yielding the probability that a model gives the set of data x. 

Suppose we want to test an hypothesis, i.e. a cosmological model set by M 
parameters = (9^, 9 2 , ... #ju). The likelihood function L(x\0) = exp[— £(x\0)] is 
often a complicated function of 6; the value corresponding to the peak of L defines 
the maximum likelihood estimator which, in the limit of large data sets, becomes the 
best unbiased estimator of the actual parameter set. Thus, the likelihood can be 
Taylor expanded to second order (the first non-vanishing term) around 6, being so 
approximated with a multivariate Gaussian distribution 



L(x\6) oc exp 



-Ae T c(e) 1 A9 

2 v ; 



here 



c(oy 



d 2 £(x\0) 



dO n dO 







(20) 



(21) 



9=9 



a/3 



the average (...) 



is a positive semi-definite non-singular matrix, dubbed covariance matrix of the 9 a . We 
remind that Equation (120]) holds just in a sufficiently small neighborhood around the 
maximum 0. In turn, the Fisher information matrix reads 

„ / d£(x\0)d£(x\0) \ _/ d 2 £(x\0) \ 

' \ d9 a d9 p J Q=e \ d9 a d9 f} / e=e ' 1 ] 

J £(x\6) . . . d N x is taken over all possible data realizations, given the 
model parameters. The feature making Fisher's formalism so significant is the Cramer- 
Rao theorem. It states that the parameter variance about any unbiased estimator value 
owns a lower bound: A9 a > a/ (F -1 )^, if the other parameters are estimated from the 
data as well, A9 a > 1/ \f~F aa , if all the other parameters are known. Therefore, the Fisher 
information matrix components are the expectation values of C _1 (0). Accordingly, the 
inverse of the Fisher matrix is an estimate of the covariance matrix of the parameters 
C(0) w F" 1 . 

A convenient way to re-write the Fisher matrix is computing the derivatives of the 
likelihood function using the following chain rule [38J: 



d£{x\0) 



39 n 



E 



d£{x\0) 



dxf 



dx £ 



X=X{0) 



e=e t 

Thus, the Fisher matrix ff22l can be expressed as: 

dxpdxpi I d 2 £(x\0) 



d9 n 



9=0 



E 

w 



d9 a 89 p 

dxo _ . 



dxedx£i 
dxoi 



X=X(9) 



(23) 
(24) 
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where and C«/, respectively, are the Fisher and the covariance matrix for the 
observables x. The region in the M-dimensional space of the parameters, defined by 
Q(9, 9) = AOaFapAOp = K 2 , is a hyper-ellipsoid of constant probability density for the 
function (|20|) . Marginalizing over the other parameters, one can project this ellipsoid 
in the two-parameter subspace, yielding a two-dimensional ellipse. The analytical 
expression for the projected ellipse for the two parameters 9 a and 9p is given by |39j : 



A9 n A9 



)aj3 




A X 2 (N = 2,a) (27) 

This can be interpreted as an estimate of the confidence region within a given confidence 
level a for the two parameters 9 a and 9 p. 



3. Forecasts for Future Experiments 

We present here the results of the Fisher analysis of future experiments, considering 
both WL and CMB measurements. For definiteness we assume a fiducial WL survey 
with characteristics similar to those of the recently proposed DUNE project [ID]. We 
assume a redshift distribution of the form (|T5|) with A = 2, B = 1.5 and z ~ z m /1.412, 
corresponding to a median redshift of the survey z m = 0.9 (see [27J), and a mean 
surface density of galaxies n g = 35arcmin~ 2 . The full survey, covering half of the sky 
(fsky = 0.5), is divided into N = 5 redshift bins, with p(z p h\z) given by equation ( fT4l) 
and a z {z) = 0.05(1 + z). 

We consider lensing multipoles up to £ max = 20000, since we find that results do 
not depend significantly on larger £. However, one should bear in mind that when 
considering scales £ ^> 1000 there could be some non-linear and baryonic effects on the 
matter power spectrum, and so on the WL spectrum |41j . These effects, not yet well 
understood, could be important for forecasts. However, in this work we suppose these 
effects to be negligible. 

For CMB data, we consider an ideal experiment with characteristics based on 
the 143GHz PLANCK channel: angular resolution #f w h m = 7.1' and sensitivity o"t = 
42 / uKarcmin, <Tp = 80;uKarcmin. 

The cosmological model we consider is characterized by 7 parameters with fiducial 
values: 

= { Uh = (0.045 • 0.7 2 ), uj m = (0.30 • 0.7 2 ), fi m = 0.30, n s = 1.00, 
cr 8 = 0.8, A DE = 5 • 10~ 3 GeV, (3 = 0.1 }. 

Here Q m represents the current total (CDM + baryons) matter density in units of the 
critical density; = fib h 2 and u m = Q m h 2 are the physical baryons and total matter 
densities, respectively; n s is the slope of the primordial power-law spectral index of 
density fluctuations; o"g is the rms mass fluctuation in spheres of 8ft. -1 Mpc radius while 
Ade and (3 were defined in Sec. 12.11 Let us notice that the class of DE models considered 
here reduces to ACDM for (A DE ) 4 ~ 10~ 47 GeV 4 and (3 = 0. Moreover, the fiducial values 
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Table 1. Estimated errors on model parameters. 

CMB WL WL+CMB 

SUGRA SUGRA SUGRA SUGRA SUGRA SUGRA 





= 0. 


g = o.i 


= 0. 


= 0.1 


= 0. 


= 0.1 


100*u b 


0.016 


0.019 


0.5 


0.9 


0.011 


0.012 




0.002 


0.006 


0.016 


0.03 


0.0004 


0.0005 




0.05 


0.12 


0.002 


0.0014 


0.0011 


0.0014 


n s 


0.004 


0.005 


0.012 


0.018 


0.0014 


0.0021 


&s 


0.07 


0.13 


0.0026 


0.0029 


0.0017 


0.0016 


A 


7.2 


9.5 


0.89 


1.1 


0.28 


0.28 


P 




0.04 




0.018 




0.0016 


T 


0.005 


0.006 











of DE parameters A^e = 5 ■ 10 _3 GeV and (3 = 0.1 are chosen in order to reproduce at 
z = an effective equation of state which mimics the case of ACDM model, w = —0.95. 
Finally, when dealing with CMB data, we also need to fix the value of the optical depth 
to reionization, r = 0.10. 

We compute the CMB anisotropies (temperature and polarisation) power spectra 
and the transfer functions, used to calculate linear matter power spectrum, using a 
modified version of CAMB [12]. To evaluate the non-linear matter power spectrum, 
Pnli we employ the prescription by Smith et al. [33]. This is only tested for model with 
a cosmological constant; as we are concerned here with Fisher matrix estimates assume 
that the results of [13] can be extended to coupled models simply by taking into account 
the non-standard scaling of p c (eq. [3]). Numerical derivatives were evaluated considering 
a 5% stepsize, except for A DE , where we adopted a 5% stepsize on A = Logi (A DE /GeV). 

3.1. CMB measurements 

Table [1] lists the estimated errors on the various parameters considered. For each data 
set, we compare forecasts for the target model with results for a SUGRA model with the 
same values of the relevant parameters. The table clearly shows that a PLANCK-like 
experiment is able to provide a measurement of a direct DE-DM interaction at 68% 
confidence level, even for moderate values of the coupling strength (3. However, we 
expect that at 90% confidence level data will still be compatible with f3 = 0. 

In any case, allowing for a direct interaction strongly degrades the experimental 
sensitivity on the parameters characterizing the matter density and the normalization 
of the primordial spectrum of density fluctuations. Errors on these quantities increase 
by a significant amount. 

Figure [2] shows the joint 68% confidence regions between f3 and each of the other 
parameters, except for r, considering only CMB data. In each plot we marginalized over 
the parameters not shown. (3 is strongly correlated with most parameters considered 
here, with the exception of uo^ (and r), thus introducing additional degeneracies in actual 
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Figure 2. Forecasts of joint 1-er confidence regions on the coupling parameter (3 = O.f 
and selected parameters, for a PLANCK-like experiments, after full marginalization 
over the remaining parameters. 



data analysis. 

A detailed characterization of these degeneracies would require a different approach 
than that followed here (e.g. Monte Carlo Markov Chains simulations). We just point 
out that they can be understood recalling that the heights of the acoustic peaks of CMB 
spectra are sensitive to the total matter density and to baryon/dark matter ratio at last 
scattering. In coupled models, these quantities are not univocally determined by their 
present day value, but strongly depend on (3. In addition, the total growth between the 
last scattering epoch and today is strongly sensitive to (3, resulting in a clear degeneracy 
between j3 and cr 8 . 

Let us notice that, as we are concerned with a combination of CMB and WL 
data, the parameter set is not optimized for CMB experiments. Using a different 
parametrisation can alter error estimates and/or the degeneracies between the various 
parameters. In particular, CMB data are better described in terms of the angle 
subtended by the acoustic horizon at recombination, 9, and the amplitude of the 
primordial spectrum of density fluctuations A s , rather than in terms of Q m and as- This 
results in large errors on the latter parameters, when CMB data alone are considered. 
Adopting a set of parameters better suited to the analysis of CMB data results in slightly 
lower error estimates overall, but the effects of coupling are largely unchanged. 

3.2. Weak Lensing 

Figure [3] is analogous to figure [2] for our target weak lensing survey. Together with 
Table [U these results show the great potential of WL surveys in constraining interacting 
DE models. Marginalized errors on f3 and A are of the order of o~{(3) ~ 0.02 and 
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Log 10 (A DE /GeV) 



Figure 3. Forecasts of joint 1-a confidence regions on the coupling parameter (3 = 0.1 
and selected parameters, for a DUNE-like experiments, after full marginalization over 
the remaining parameters. Notice the change of scales with respect to figure [2] 



0.125 - 



0.1 - 



0.075 - 



0.05 - 



0.025 - 








-4 -2 
L°g,„ (A DE /GeV) 



Figure 4. Comparison between the 1— a confidence regions of a coupled SUGRA model 
with (3 = 0.1 (red) and a non-coupled SUGRA model with (3 = (blue) for a WL 
experiment. It is clearly possible to distinguish between the two models. 

<t(A) ~ 1; these figures represent a factor of 2, or more, improvement over Planck 
estimates. WL data alone can clearly distinguish the target model from a non-coupled 
model or a cosmological constant even at the 3cr-level, viceversa assuming a reference 
SUGRA model with /3 — 0, we can expect to put an upper limit (3 < 0.03, at the same 
confidence level (see figure SJ) . 

As expected, WL surveys perform significantly better than CMB experiments also 
with respect to parameters specifying the current matter density, Q m and o~&. Moreover, 
constraints on these parameters are not significantly affected by the coupling degrees 
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-4 -2 0-4 -2 

Log 10 (A DE /GeV) Log 10 (A DE /GeV) 

Figure 5. Joint 1-tr confidence regions on /3 and Ade after marginalization over the 
remaining parameters, for different number of galaxies n g . On the left SUGRA coupled 
model with j3 = 0.1, on the right SUGRA coupled model with (3 = and j3 derivatives 
calculated only on one side, for positive values of the parameter. 




1 1.5 2 2.5 3 3.5 4 4.5 1 1.5 2 2.5 3 3.5 4 4.5 

L°6,„(l) Log,„(l) 



Figure 6. Inverse error as a function of the maximum multipole for the WL survey, 
for different numbers of redshift bins. Left panes show results for /3, right panel refers 
to A. 

of freedom. Errors on the remaining parameters, instead, increase by a factor of ~ 2. 
Finally, we consider a combination of CMB and WL data. CMB and WL probe very 
different epoch of the Universe and are sensitive to different combination of cosmological 
parameters. Considering both CMB and WL data allows to constrain the DE parameters 
with a few percent accuracy, and significantly reduces the degeneracies introduced by 
DE coupling. In this case, the errors on the cosmological parameters are very similar in 
both models considered, with the exception of n s . 

Next we considered how our results depend on the characteristic assumed for 
the target survey. In figure [5] the impact of mean surface density of galaxies on the 
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determination of f3 and A. With n g = 25 arcmin -2 constraints on (3 degrades by 
~ 50%, while n g = 50 arcmin -2 gives only a marginal improvement on expected errors; 
constraints on A are similarly affected. However, even in the worst case considered here, 
next generations WL survey will provide an improvement over the information that we 
are likely to obtain from PLANCK data. 

Lastly, we consider the dependence of our results on the number of bins and the 
multipoles range considered. In figure [6] we plot the inverse of the expected variance 
on (3 and A as a function of the maximum multipole considered in the analysis and for 
different number of bins. With 3 redshift bins, the precision on both parameters depends 
mostly on multipole up to a few thousands; smaller scales do not provide a significant 
contribution. Dividing the survey in 5 bins strongly improves the constraints on both 
parameters and allows to exploit information from multipoles up to ~ 10000. For a 
DUNE-like survey, a further increase in the number of bins does not lead to significant 
improvements on the constrains on coupled models parameters. 

It must be outlined that these results assume that the theoretical framework used 
to predict the matter power spectrum on intermediate (1/iMpc -1 < k < 20/iMpc -1 ) and 
small (k > 20/iMpc -1 ) scales can accurately account for the effects of baryons on non- 
linear structures. In general, the fitting formulas used to predict the non-linear power 
spectrum are calibrated using dissipationless N-body simulations and, therefore, do not 
properly describe baryonic structures. While baryons make up ~ 15 — 20% of the matter 
in the Universe and on large scales are expected to trace the DM field, their distribution 
inside halos is significantly different from DM. In turn, this alters the shape of the non- 
linear power spectrum on the corresponding scales, and the possibility of extracting 
precision constraints from Pnl(^) hinges on our capability of accurately modeling baryon 
physics [41J. However, simulations do not yet have the accuracy required for precision 
constraints and the problem is even more serious for the coupled models considered in 
this work. Modelling non linear stages though spherical growth, Mainini |44j shew that 
baryons and DM will be however differently distributed, even independently of the onset 
of gas dynamics. N-body simulations of cDE models were performed [15] . by using a 
Ratra-Peebles [46J potential; hydro simulations, instead, were never produced. Should 
accurate prediction be still unavailable for the analysis of a DUNE-like experiment, a 
more conservative cutoff of I ~ 1000 would be required. Figure [6] shows that in this case 
the expected errors on A and (3 would increase by a factor ~ 2. 

4. Discussion 

All previous analysis shows that, even if we admit quite a little DM-DE coupling, 
we open a Pandora's box, leading to a severe degradation in our capacity to deduce 
cosmological parameters from a given set of measures. 

As a matter of fact, coupling destroys our trust that the period between the 
recombination and the start of DE relevance is under control. If coupling is absent, 
during such period SCDM is a fair approximation. Let us then remind what happens 
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Figure 7. Scale dependence of comoving distances in coupled or uncoupled SUGRA 
cosmologies, compared with ACDM. The case (3 = 0.5, corresponding to a rather strong 
DM-DE coupling is also shown. 



to the growth factor, as soon as coupling is onset: Figure ED shows that: (i) deviations 
from SCDM are significant already when a ~ 0.02; (ii) they are then different for DM 
and baryons; (iii) they work in the opposite direction, in respect to the effects of a DE 
components. 

As far as the growth factor is concerned, a tiny coupling is able to overwhelm a 
huge DE amount, with compensation occurring for a ~ 0.3-0.4 , however keeping always 
g(a) at values greater by ~ 10-15 % . Altogether, growth is faster in coupled models. 
Hence, if we do not include the information that coupling is zero in the fit, we can find 
an agreement between data and a wider range of DE amounts. 

Similar points can be made for the comoving radial distance R(z) = r(0, z) (see 
eq. [9]). In Figure [7] we compare comoving distances for ACDM with various cosmologies. 
The Figure shows that dDE, in the absence of other parameter shifts, sets the Last 
Scattering Band (LSB) closer to the observer. Once again, a mild coupling acts in 
the opposite direction and tends to re-set the LSB at the distance it had in ACDM. 
In the Figure we consider the behavior of distances also for a rather strong coupling, 
(3 = 0.5. Then the distance behavior is different in the period when DE density can be 
neglected, in respect to the epoch when DE and DM have similar densities. The key 
point, however, is that the LSB becomes then farther from the observer. When fitting 
CMB data to such models, in order to compensate such effect, the value of H Q tends to 
be increased. Strong coupling therefore yields a large Hubble parameter estimate. 

Figures [8] finally show the scale dependence of the density parameters in the different 
models. Once again, when DE is mildly coupled, a behavior more similar to ACDM is 
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recovered. On the contrary, when considering a greater coupling strength, we see that 
DE and DM keep similar densities up to a fairly large redshift. This was indeed the 
initial motivation of cDE cosmologies. 

Altogether, these Figures indicate that adding a small coupling reduces the effects 
of the very passage from ACDM to dDE; owing to the excellent fit that ACDM 
cosmologies have with data, this tells us that only highly refined CMB data will be 
able to test the possibility that a mild DM-DE coupling exists. 

On the contrary, a stronger coupling, although easing the coincidence problem, 
displaces several observables in a unacceptable way. 

5. Summary and Conclusions 

Future WL surveys will certainly put more stringent constraints on cosmological 
parameters and will be crucial to break quite a few degeneracies between parameters. 

Within this context, in this paper we focus on coupled DE models with a twofold 
aim. Detecting a signal of DM-DE coupling would be certainly decisive to fix the 
nature of the dark components. Henceforth determining the level of sensitivity needed 
to appreciate such an effect is crucial in setting the appeal of forthcoming projects. 
There is however a complementary aspect which deserves much attention. In order to 
convert raw data into physical information, a set of parameters, spanning a variety of 
models, is to be fixed; a bias on parameter selection, however, can lead to an optimistic 
estimate of the confidence level for the best fitting model, far from reality. 

In this paper, we focused on this kind of danger, when we open the option of 
DM-DE coupling. Even if such coupling is absent or quite weak, we showed that just 
considering its possibility may widen the error bars for a number of parameters, also 
apparently unrelated to the coupling itself. In a sense, when a new degree of freedom 
is opened, such an effect is natural and expected. Coupling, however, has really a 
major impact, affecting different parameters for the different observables considered 
here; moreover, its impact is drastically reduced when we work out parameter values by 
using simultaneously both observables. 

Before outlining our main conclusions it is however worth reminding a technical 
limit we had to face. While CMB predictions depend on linear spectra, the WL spectrum 
is limited to its non-linear shape. For the purpose of the present analysis, we assumed 
that prescription for ACDM models [33] can be trivially extrapolated to coupled models, 
so enabling us to estimate the non linear spectrum once the linear spectrum is known. 
Let us however remark that the shift was estimated from quite a wide set of ACDM 
simulations. Simulations of cosmological models with state parameters w^l, although 
performed by several authors, are still not so extensively studied as ACDM; let alone 
coupled DE simulations: in this case the only available simulations are due to [35] and 
deal with a potential V((f>) different from SUGRA. 

However, differences between prescriptions for ACDM and dDE are small [37J and it 
seems however clear that model differences can most affect the rate of evolution of halo 
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coupled SUGRA models are considered. For the sake of comparison, in the lower panel 
we also show the effects of choosing a stronger coupling. 



concentration, slightly shifting the scale where non-linearity effects become significant. 
The use of more precise prescriptions can therefore only cause minor variation on the 
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Figure 9. Forecasts of joint 1-a confidence regions on the coupling parameter = 0.1 
and selected parameters, for a combination of a PLANCK-like and a DUNE-like 
experiment, after full marginalization over the remaining parameters. 

estimated errors and, although welcome, such simulations are not expected to interfere 
substantially with our conclusions. 

Our estimates were based on assuming that a photometric survey is available, with 
o~z(z) = 0.05(1 + 2;) and a median redshift z m = 0.9, covering half of the sky (f s k y = 0.5). 
These features are similar to the recently proposed DUNE experiment [ID]. We also 
compared and combined results from WL with the constraints expected for an ideal 
Planck-like experiment. The basic results of our calculations are quoted in table Q] 
which is one of the main results of this work. 

A first set of conclusions concerns an ideal CMB experiment considered by itself. 
In this case, introducing coupling degrees of freedom is crucial for the error estimates 
on some of parameters, also apparently unrelated to coupling. In particular, while the 
error on keeps ~1%, the errors on u m and <Tg increase from 1.3% to 3.9% and from 
9% to 16%, respectively. 

From a physical point of view, the option opened by coupling is that expansion 
rate and fluctuation growth, from the last scattering band to the observer's site, is non- 
standard {e.g., the proportionality law p m oc a~ 3 could be mildly violated). Although 
CMB data themselves set stringent limits on such deviations, this widens the volume 
of the parameter space consistent with a given data set; in particular, it increases 
the likelihood of values of Ho that would otherwise be negligible and, because of the 
intercorrelation amongst parameters, this reflects immediately on uo m and Q m estimates. 

Similar effects occur in WL experiments, although involving different parameters. 
Such experiments are a direct test of fl m , whose estimated error is reduced by a factor 
> 30 in respect to a CMB experiment. When the coupling option is opened, the error 
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on Q m does not increase; on the contrary, it becomes easier to attribute raw data 
uncertainties to other parameters and, although marginally, the error on Q m becomes 
somehow smaller. 

A completely new situation occurs if both CMB and WL measurements are 
simultaneously used. In this case, the opening of the coupling option causes just a 
marginal increase of the errors on most parameters. This is a clear indication of the 
complementarity of the CMB and WL measurements, as described in [H] and one of our 
conclusion is that the combination of these observables, besides of providing parameter 
values independent of /3, can set a (nearly-) final word on the coupling option. 

As a matter of fact, by comparing figure El figure [3] and figure [HJ we see that, when 
joining CMB and WL results, the degeneracies, between (3 and o~ 8 , as well as between f3 
and fl m , disappear. Breaking degeneracies is the main aim when different observables 
are simultaneously considered. We see that, from this point of view, the efficiency of 
using both CMB and WL measures can be hardly overestimated. 

Let us then focus on the case of the spectral index of scalar fluctuations. When 
both CMB and WL data are used to constrain n s , a sharp reduction of errors occurs. 
No surprise that CMB data, by themselves, reflecting the state of the universe before 
the onset of non linear processes, were more efficient to constrain n s , with or without 
coupling, than WL measures. Joining together the two observables, we then see errors 
to decrease from 0.4% and 1.2% down to 0.16%, in the uncoupled case; in the presence 
of coupling we have a similar behavior, with errors passing from 0.5% and 1.8% down 
to 0.25%. The error level achieved, in both cases, is exceptional, even for precision 
cosmology, and clearly suggests to relieve the constraint of a single n s value, so inspecting 
its possible scale dependence, with realistic possibilities to find a direct insight into the 
nature of the inflationary potential. 

In conclusion, future WL surveys could really allow a significant step forward in the 
comprehension of the dark cosmic side; we can adfirm that, when they will be available, 
the endeavour to put the genie back inside the lamp will approach a full success. 
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Appendix A. Convergence power spectrum covariance 

In order to determine the convergence power spectrum covariance, one can introduce the 
so called "flat-sky" approximation and treat the sky as flat, replacing spherical harmonic 
sums with Fourier transforms (FT). Of course, this approximation is acceptable just for 
small angular scales We also consider the tomographic case and use Greek letters as 
superscripts to denote quantities belonging to different redshift bins. The FT of the 
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convergence field can be denned as: 

K a (1) = J d 2 K a {6) exp {-W ■ 1), 

while the convergence power spectrum and trispectrum are, respectively: 
(K a (h) K P (h)) = (2n) 2 6 D (l 1 + h)P^, 

{K a . . . k 5 (h)) c = (2vr) 2 5 D (h + ... + 10 T a ^ s (h, l 2 , l 3 , 1 4 ) ; 
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(A.l) 



(A.2) 
(A.3) 



here 5d is the Dirac function. The value of the lensing power spectrum in correspondence 
of a multipole I can be estimated as the mean over a multipole bin of width Al, centered 
at I: 



' l ~ A 



d 2 lx 

Ai 



K a (10 K? (-li 



(A.4) 



where A t = J s d 2 l = 2vr/ Al is the area of the shell of width Al corresponding to I, while 
A = 4n f sky is the area of the survey. Quite in the same way, for the trispectrum we 
have: 



(A.5) 



Let us then consider the following expression: 



{vfvf) 



A 2 



d 2 l! 
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The 4-point function at the r.h.s. can be decomposed in its connected parts: 

(«« (10 J (-10 (10 k 5 (-10) = (k« (10 ^ (-10 ^ 7 (la) (-Ia)>c + 

+(k«(10^(-10)(k 7 (1 2 )k 5 (-1 2 )) + 

+<k« (10 (1 2 ))(^ (-10 (-I2)) + (k q (10 (-h)}(^ (h) k? (-10). (A.7) 

Replacing their expression in eq. (1A.6I) . one can easily recognize the contribution of the 
trispectrum, using eq. (1A.3I) : 



(VfV?, s ) = j- 2 f^f ^(2tt) 2 ^(0)T^(1 1 ,-1 1 ,1 2 ,-10 + (^)(^) 
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( K a (10 ^ (i 2 ))(^(-i0^(-i 2 )) + 
<«■ (io^(-iO)(^(i 2 )^(-iO), 



where (27r) 2 5d (0) = A. Owing to the definition of covariance, 



Cov 



1 I 1 'V 



{vfvf)-{vr){vT), 



(A.8) 
(A.9) 

(A.10) 



and using (1A.2j) . one can then argue that the integrals in (1A.8j) make sense only if they 
correspond to the same /-bin; the same can be claimed for ( 1A.9I) . This property can 
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be described introducing a Kronecker delta function Sur. Thus, the expression (IA. 10h 
becomes: 
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4:71 f sky 2lAlf sky 

In the last line, we have supposed the lensing power spectrum to be smooth enough to 
treat it clS cl constant within each bin width. 



(A.12) 
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